Late-Time Evolution of Charged Gravitational Collapse and 
Decay of Charged Scalar Hair - III. Nonlinear Analysis 



Shahar Hod and Tsvi Piran 

The Racah Institute for Physics, The Hebrew University, Jerusalem 91904, Israel 

(February 7, 2008) 

Abstract 

We study the nonlinear gravitational collapse of a charged massless scalar- 
field. We confirm the existence of oscillatory inverse power-law tails along 
future timelike infinity, future null infinity and along the future outer-horizon. 
The nonlinear dumping exponents are in excellent agreement with the ana- 
lytically predicted ones. Our results prove the analytic conjecture according 
to which a charged hair decays slower than a neutral one and also suggest 
the occurrence of mass-inflation along the Cauchy horizon of a dynamically 
formed charged black-hole. 

I. INTRODUCTION 

Linearized perturbation analysis have revealed important dynamical features of the grav- 
itational collapse of massless neutral fields. In particular, according to the linearized per- 
turbation theory there are two major features which characterize the late stages of the 
evolution: quasinormal ringing and inverse power-law tails (which follows the quasinormal 
ringing). These late-time tails are relevant to two major aspects of black- hole physics: The 
no-hair theorem of Wheeler and the mass-inflation scenario [Q . The mechanism responsible 
for the development of these neutral inverse power-law tails was first studied by Price [0. 
This work was further extended by Gundlach, Price and Pullin Q. These authors have 
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also shown that the nonlinear dumping exponents, describing the fall-off of a massless neu- 
tral scalar-field at late-times, are in good agreement with the prediction of the linearized 
perturbation analysis ||. 

The late-time evolution of a charged massless scalar-field and the physical mechanism for 
the radiation of a charged hair were first studied analytically in hereafter referred to as 
papers I and II, respectively. The main result presented there is the existence of oscillatory 
inverse power-law tails along the asymptotic regions of future timelike infinity, future null 
infinity and along the future outer-horizon with smaller (compared to the neutral case) 
dumping exponents. These dumping exponents depend on the black-hole's (or the star's) 
charge. In this paper we study numerically the fully nonlinear gravitational collapse of a 
massless charged scalar-field. We confirm that the late-time behaviour of the fully nonlinear 
evolution is in excellent agreement with the analytical predictions of the linearized analysis. 

The plan of the paper is as follows. In Sec. we give a short review of the linearized 
analytical results of papers I and II, on which we base our expectations for the fully nonlinear 
case. In Sec. [ITJ we describe our physical system, namely the coupled Einstein-Maxwell- 
charged scalar equations. In Sec. [TV] we study the late-time evolution of the nonlinear 
spherical charged gravitational collapse. In Sec. |V| we study the late-time evolution of non- 
spherical perturbations on a fixed Reissner- Nordstrom background and on a time- dependent 
background. In all cases we compare the nonlinear results with the predictions of the 



linearized analytical theory. We conclude in Sec. |VI] with a brief summary of our results 
and their physical implications. 



II. REVIEW OF LINEARIZED ANALYTICAL RESULTS 

We study numerically the late-time behaviour of the fully nonlinear gravitational collapse 
of a massless charged scalar-field. Our expectations are based on the linearized analytical 
results of papers I and II. Due to the smallness of the field's amplitude at late times, we 
expect these results to hold even in the fully nonlinear case. In other wards, we expect 
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that the late-time field may be regarded as a small perturbation. Thus, our expectations 
include the existence of inverse power-law tails along the three asymptotic regions: timelike 
infinity i + , future null infinity scri + and along the black- hole outer-horizon H + (where the 
power-law is multiplied by a periodic term). Quantitatively, in paper II it was found that 
the late-time behaviour of linearized charged-scalar perturbations (on a Reissner- Nordstrom 
background) is dominated by power-law tails of the form 

*P l m = Ml,eQ)y^r^ , (1) 

at timelike infinity i + , 

ii l m = A ni (l,eQ)v; ie %;V +1 - ie W , (2) 
at future null infinity scri + and 

*Pl = A eh (l,eQ)e^ y v;^ , (3) 
along the black- hole outer- horizon H + , where the charged scalar-field is given by 

<P = e^ t Y,^m(t,r)Yr(9^)/r, (4) 

l,m 

where <£> is merely a gauge constant of the electromagnetic potential A t , i.e. its value at 
infinity (A t is given by A t = $ — ^). The constant (3(1, eQ) is given by 

^_ -l + V / (2/ + l) 2 -4(eg)^ (s) 

The tortoise radial coordinate y is defined by dy = dr/X 2 , where A 2 = 1 — ^ + Here 
v e = t+y and u e = t— y are the ingoing and outgoing Eddington-Finkelstein null coordinates, 
respectively. The expressions for the coefficients A(l, eQ) are given in paper II. 

III. THE EINSTEIN-MAXWELL-CHARGED SCALAR EQUATIONS 



We consider a spherically-symmetric self-gravitating charged scalar-field 0. The system 
is described by the coupled Einstein-Maxwell-charged scalar equations. This physical system 
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was already described in a previous paper j7|. Here we give only a short description of the 
system and the final form of the equations studied. We express the metric of a spherically 
symmetric spacetime in the form P,p| 



ds 2 = —g(u, r)g(u, r)du l — 2g(u, r)dudr + r 2 dfl 2 



(6) 



in which u is a retarded time null coordinate and the radial coordinate r is a geometric 
quantity which directly measures surface area. The coordinates have been normalized so 
that u is the proper time on the r = central world line. We use the auxiliary field <fi defined 
by 



r Jo 



dr 



(7) 



in terms of which the Einstein equations yield 



g(u, r) = exp 



Art 



dr 



(8) 



and 



The electromagnetic potential A u is given by the Maxwell equations 

r Q 



(9) 



A, 



o r z 



-g dr 



(10) 



where the charge Q(u,r) within a sphere of radius r, at a retarded time u is given by 



Q(u,r) = Airie / r(0*0 — (jxf)*)dr . 
Jo 

The mass M(u, r) within a sphere of radius r, at a retarded time u is given by 



(11) 



M(u,r) 



2 7r £(0_0)(J*_^) + ±5L 

g ) ) ■ - 



1_Q 2 

2 r' 



dr+ 2V 



(12) 



Finally, the wave-equation for the charged scalar-field takes the form of a pair of coupled 
differential equations 



dcj) 
du 



1 



(g-g)(<j>-4>) 



Q 2 



(4>-4>)g 



ieQ 



g<p — ie(f)A. 



(13) 



2r 



2r 3 



2r 



M 5 



and 




1 



(14) 



We solve this system of equations numerically. A detailed description of our algorithm, 



In order to compare our nonlinear results with the linearized analytical results of papers 
I and II, we use another spacetime coordinate, namely the Bondi time ts = t(u, oo), which 
is the retarded time coordinate that agrees with time at infinity. The proper time t(u, r) 



The comparison with the analytical results also requires the usage of the ingoing and outgo- 
ing Eddington-Finkelstein null coordinates v e and u e . In the asymptotic region r 3> Mbh-, 
v e is linear with r (along a« = const ray). In particular, v e ~ 2y ~ 2r along the u = u e = 
ray. In a similar way u e is linear with t along a v = const ray, namely u e = 2t — v e (one may 
also use the asymptotic relation u e ~ v e — 2r along a v = const ray in order to evaluate u e ). 
The relation between A u and A t implies that A t {r + ) = 0, i.e. $ = 



In this section we present our results for the nonlinear spherically symmetric gravitational 
collapse of the self-gravitating massless charged scalar-field. We have focused our attention 
on the behaviour of the charged field <fi along the three asymptotic regions: timelike infinity 
i + , null infinity scri + and the black- hole outer horizon H + . The late-time evolution of a 
charged scalar-field is independent of the form of the initial data. The numerical results 
presented here have an initial profile of the form 



numerical methods, discretization and error analysis are given in Ref. 0]. 



along an r = const trajectory is given by [|j] 




(15) 



IV. NONLINEAR SPHERICAL COLLAPSE 
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4>(u = 0, r) = Ar 2 exp j— [(r — To)/ "] 2 } > (16) 

where ro = 2.5, cr = 0.5 for the real part of the complex field and r$ = 3.0, a = 0.5 for 
the imaginary part. Figure 1 displays the time evolution of 4>r, the real part of the charged 
scalar-field, along these asymptotic regions for a supercritical evolution (A = 0.005, e = 0.85) 
in which a black-hole forms. The mass and charge of the formed black-hole are Mbh — 0.503 
and Qbh — —0.420, respectively. The top panel shows the behaviour of the field at a con- 
stant radius (here r = 10) as a function of the Bondi time. The middle panel shows the 
behaviour of the field along null infinity (approximated by the null surface v = v max , where 
v max is the largest value of v on the grid) as a function of u e . The bottom panel shows the 
behaviour of the field along the black-hole outer horizon H + (approximated by the null sur- 
face u = u max , where u max is the largest value of u on the grid) as a function of v e . Initially, 
the evolution is dominated by the prompt contribution and by the quasinormal ringing. 
However, at late-times a definite oscillatory power-law fall off is manifest. The nonlinear 
power-law exponents (determined from the maximas of these oscillations) are —1.78 at time- 
like infinity, —0.86 along null infinity and —1.96 along the black-hole outer horizon. These 
values should be compared with the analytically predicted values of —1.70,-0.85,-1.70, 
[see Eqs. ([[]) — (||)], respectively. The theoretical values of the oscillation frequencies of the 
charged scalar field are at timelike infinity and along the black-hole outer horizon 

and a l° n g nu ll infinity [see Eq. @ and the relation $ = ^-]. These values agree to 

within 2% with the numerical estimates. 

While the dumping exponents and the decay-rate of neutral perturbations are indepen- 
dent of the spacetime parameters (M and Q), the analysis of paper II predicts that the 
charged dumping exponents depend on the black-hole's (or the star's) charge (namely, on 
the dimensionless quantity |eQ|). In order to test this linearized analytical prediction, we 
have studied the dependence of the charged dumping exponents at timelike infinity on the 
parameter e. The initial form of the field is given by (16) and the amplitude is set to 
A = 0.005. The nonlinear results are summarized in table |. We find a good agreement 
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between the numerically measured exponents and the linearized predicted ones. 

The analysis of paper I predicts the existence of charged power-law tails even in subcrit- 
ical evolution when there is no collapse to a black-hole (and all we have are imploding and 
exploding shells). This phenomena is related to the fact that the late-time evolution of the 
field is dominated by the backscattering from asymptotically far regions, and it does not 
depend on the small-r details of the spacetime (this is also the situation for neutral massless 
perturbations ||). Since the charge of the spacetime is a dynamical quantity it is not clear 
which value of Q should be taken in calculating the value of the dumping exponent. However, 
since the charge of the spacetime falls to zero asymptotically we expect the effective value 
of \eQ\ to be much smaller than unity. In this case the dumping exponents are expected to 
be 21 + 2 at timelike infinity and I + 1 along null infinity (with small corrections of order 
0[(eQ) 2 }). Figure 2 displays the time evolution of 4>r for a subcritical evolution (A=0.002, 
e=3). Shown are the behaviour of the field at a constant radius (here r = 10 ) as a function 
of Bondi time and along null infinity as a function of u e . The nonlinear power-law exponents 
are —2.00 at timelike infinity and —1.00 along null infinity. These values are exactly equal 
to those predicted in paper I. 



V. NON-SPHERICAL PERTURBATIONS OF SPHERICAL COLLAPSE 



In order to study the dependence of the late-time dumping-exponents on the multipole 
index I we have performed two series of numerical investigations. First, we have integrated 
the linearized charged scalar-field equation 



Q 

tp jt t + 2ie—il), t - ip m + Vip = 



(17) 



where 



V = V r 



M,Q,l 



I (I + 1) 2M 2Q 2 



— e 



2 Q 2 



(18) 



on a fixed Reissner-Nordstrom background. It is straightforward to integrate Eq. ( |T?D using 
the method described in ||. The late-time evolution of a charged scalar- field is independent 



of the form of the initial data. The results presented here are of a Gaussian pulse on u = 

rj)(u = 0,v) = Aexp[- [(v - v )/a] 2 } , (19) 

where the amplitude A is physically irrelevant due to the linearity of Eq. (|T7|). It should be 
noted that the evolution equation ( |TTD is invariant under the rescaling 

r — > ar , t — > at , M — > aM , Q — > aQ , e — > e/a , (20) 

where a is some positive constant. The black- hole mass and charge are set equal to Mbh = 
0.5 and Qbh = 0.45, respectively. We have chosen e = 0.01 , v o = 100 and a = 20 (for 
both the real and the imaginary parts of the charged field). The numerical results for the 
/=0,1 and 2 modes are shown in Fig. [| (from top to bottom, respectively). This figure 
demonstrates the dependence of the late-time tails at i + (here y=400) on the multipole 
index /. A definite power-law fall off is manifest at late-times. The numerical values of the 
power-law exponents, describing the fall-off of the field at late times are —1.97, —3.94 and 
— 5.75 for Z=0,1 and 2, respectively. These values are to be compared with the analytically 
predicted values of —2.0, —4.0 and —6.0, respectively. 

In a second series of numerical investigation, we have studied the evolution of a sec- 
ond perturbative charged scalar-field £ evolved on a time- dependent spacetime. This time- 
dependent spacetime is determined by the background solution <ft, while we ignore the contri- 
bution of the field £ to the energy-momentum tensor (this is analogous to the case studied in 
0j for a neutral field). Resolving the field into spherical harmonics £ = £L(^> r )^ m (^) 

l,m 

(and using the spherical symmetry of the time-dependent background) one obtains a wave- 
equation for each multiple moment of the field £ 

% = Yr {9 - ^ - - ~ ej9 ~ 1?^™ ~ i6 ^ Au ~ Yr l{l + 1)9 ^ ' (21) 

[together with Eq. flU])], where 

i l m = -fLdr. (22) 
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The quantities g, g and Q are still determined by the background solution for the field <j). For 
the background field <fi we choose the same initial form as in Sec. [TV]. The most interesting 
time-dependent backgrounds are the non-collapsing ones. Thus, we take the subcritical 
initial conditions ^4=0. 002 and e=3 of Sec. [IV]. For the perturbation field £ we choose initial 
data of the form Eq. ( |T6| ) with r = 4.0 and a = 0.5 for both the real and the imaginary 
parts of the field (again, the amplitude of the perturbation field is physically irrelevant). 
Figure 4 displays the time evolution of the perturbative charged scalar field £ (its real part) 
at a constant radius (here r = 10) as a function of Bondi time. It is interesting that even 
for time- dependent spacetimes the late-time behaviour of charged-fields is well described by 
an inverse power-law fall-off. The numerical values of the power-law exponents are —1.99 
for the /=0 mode (top panel) and —4.16 for the 1=1 mode (bottom panel). These values 
should be compared with the fixed background analytically predicted values of —2.0 and 
—4.0, respectively. 



VI. SUMMARY AND PHYSICAL IMPLICATIONS 

We have studied the nonlinear gravitational collapse of a massless charged scalar-field. 
Following the predictions of the linearized analytical theory (papers I and II) we have focused 
attention on the asymptotic late-time evolution of the charged-field. Our main results and 
their physical implications are: 

We have confirmed the existence of oscillatory inverse power-law tails in a collapsing 
spacetime along the asymptotic regions of future timelike infinity i + , future null infinity 
scri + and along the future outer-horizon H + . The nonlinear dumping exponents are in 
excellent agreement with the analytically predicted ones. In particular, we have verified the 
analytically conjectured dependence of the charged dumping exponents on the dimensionless 
quantity \eQ\. This dependence on the spacetime charge is contrasted to neutral pertur- 
bations, where the dumping exponents are fixed integers which does not depend on the 
spacetime parameters (namely, they are functions of the multipole index / only). We have 
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confirmed the existence of power-law tails even in non-collapsing spacetimes (i.e. imploding 
and exploding shells). 

We have also studied the late-time evolution of non-spherical perturbations (charged test 
fields) on a fixed Reissner-Nordstrom background and on a time-dependent background. In 
both cases we have found that the dumping exponents, describing the fall-off of the charged 
field at late-times, are in good agreement with the predictions of the linearized analytical 
theory (In particular, we have verified the functional dependence of the dumping exponents 
on the multipole index I). 

Our nonlinear results verify the analytic conjecture that a charged hair decays slower 
than a neutral one. Furthermore, the existence of oscillatory inverse power-law tails along 
the black-hole outer horizon suggests the occurrence of the well-known phenomena of mass- 
inflation |J along the Cauchy horizon of a dynamically formed charged black-hole. 

In a forthcoming paper we study the fully nonlinear gravitational collapse of a charged 
scalar-field, using a different numerical scheme which is based on double null coordinates. 
This scheme allows us to start with regular initial conditions (at approximately past null 
infinity), calculating the formation of the regular black- hole's event horizon, and continue the 
evolution all the way inside the black-hole. Thus, this numerical scheme makes it possible 
to test the mass-inflation conjecture during the gravitational collapse of a charged scalar- 
field. To our knowledge, the mass-inflation scenario has never been demonstrated explicitly 



before in collapsing situations (The numerical work of Brady and Smith |T(J begins on a 



Reissner-Nordstrom spacetime and the black- hole formation itself was not calculated there). 
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FIGURES 
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FIG. 1. Supercritical evolution of the charged field \(/)r\ along the asymptotic regions of time- 
like infinity (top panel), null infinity (middle panel) and the black- hole outer horizon (bottom 
panel). The initial data is a Gaussian distribution (multiplied by an r 2 factor). The field's am- 
plitude is A=0.005 and e=0.85. The mass and charge of the formed black-hole are Mbh = 0.503 
and Qbh = —0.420, respectively. A definite oscillatory power-law fall-off is manifest at late-times. 
The nonlinear power-law exponents are —1.78 at timelike infinity (here r = 10), —0.86 along null 
infinity and —1.96 along the black-hole outer horizon. These values are to be compared with the 
analytically predicted values of —1.70, —0.85 and —1.70, respectively. The oscillation frequency of 
the charged field <fi agrees with the analytically predicted value to within 2%. 
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FIG. 2. Subcritical evolution of the charged field \4>r\. The initial form of the field is the 
same as in Fig. |]. The field's amplitude is yl=0.002 and e=3. The field at future timelike infinity 
(r = 10) is shown as a function of ts- Along null infinity the field is shown as a function of u e . A 
definite power-law fall-off is manifest at late-times. The nonlinear power-law exponents are —2.00 
at timelike infinity and —1.00 along null infinity. These values are exactly the ones predicted by 
the linearized analytical approach. 
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FIG. 3. Evolution of the charged field \tpR(y = 400, t)\ on a fixed Reissner-Nordstrdm back- 
ground for different multipoles £=0,1 and 2 (from top to bottom, respectively). The black-hole 
mass and charge are set equal to Mbh = 0.5 and Qbh = 0.45, respectively. The field's charge is 
e = 0.01. The initial data is a Gaussian distribution with vq = 100 and a = 20 (for both the real 
and the imaginary parts of the field). A definite power-law fall-off is manifest at late-times. The 
power-law exponents are —1.97, —3.94 and —5.75 for the 1=0,1 and 2 modes, respectively. These 
values are to be compared with the analytically predicted values of —2.0, —4.0 and —6.0. 



14 



I I I I I I I L 

2 4 6 8 

ln(t B ) 



-4 




1.5 2 2.5 3 

ln(t B ) 

FIG. 4. Evolution of the charged test field \£,n(r = KM) I on a time- dependent spacetime. The 
initial data for the background field are those of Fig. ||. (non-collapsing case). The initial data 
for the test field £ is a Gaussian distribution (multiplied by an r 2 factor). A definite power-law 
fall-off is manifest at late-times. The power- law exponents are —1.99 for the Z=0 mode (top panel) 
and —4.16 for the 1=1 mode (bottom panel). These values are to be compared with the analytically 
predicted values of —2.0 and —4.0, respectively. 



15 



TABLES 

TABLE I. Dependence of the dumping exponents at timelike infinity i + on eQ 



e Qbh 2(3 + 2 Nonlinear exponents 

0.50 -0.241 -1.97 -1.99 

0.70 -0.344 -1.88 -1.93 

0.85 -0.420 -1.70 -1.78 

0.90 -0.443 -1.60 -1.65 
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